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In this paper I review the status of lattice QCD calculations of D and D s meson decay constants and of D meson 
semileptonic decay form factors. I restrict my discussion to results obtained from simulations with nj = 2 + 1 
sea quarks. 
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1. Introduction and Motivation 



Lattice QCD is the only systematically improvable 
calculational tool we have for quantitatively under- 
standing nonperturbative QCD effects. Accurate the- 
oretical calculations of nonperturbative QCD effects 
are essential for the experimental flavor physics pro- 
gram. One set of goals of the experimental program 
are accurate determinations of the CKM matrix ele- 
ments. This is illustrated for the weak decay process 
D — > Klv. The experimentally measured (differen- 
tial) decay rate can be written as 



oT 
If 



(known)|K,r/i(g 



(1) 



where /+(g 2 ) is one of the hadronic form factors which 
parameterize the hadronic matrix element for this pro- 
cess, (K\Vfj,\D). Hence, to determine \V CS \ from exper- 
imental measurements, we need a theoretical calcula- 
tion of the form factor with matching precision. An- 
other set of goals is to constrain beyond the standard 
model theories and to search for new physics signals. 
This effort complements the experiments at the high 
energy frontier. Accurate theoretical calculations are 
again essential. 

Since lattice QCD calculations are complicated and 
time consuming, a third important goal specifically of 
the charm physics experiments is to test lattice QCD 
methods. For example, we can use Eq.[T]to determine 
the form factors from experimental measurements af- 
ter taking | V cs | from other sources. These tests are im- 
portant to establish lattice methods for the B meson 
system, where the CKM matrix elements are less well 
known, and where input from lattice QCD is essen- 
tial. The leptonic and semileptonic D meson decays 
discussed in this talk are ideal for this. They are not 
expected to be sensitive to new physics, and the cor- 
responding hadronic matrix elements are straightfor- 
ward to calculate. Once established, lattice QCD to- 
gether with the experimental measurements can then 
be used to improve the determinations of the CKM 
angles V c d and V cs . All of these goals require accurate 
measurements and calculations. 



1.1. Introduction to Lattice QCD 

In lattice field theory, the space-time continuum is 
replaced by a discrete lattice. (For reviews of lat- 
tice QCD see Ref. This implies that derivatives 
are replaced by discrete differences, which in turn in- 
troduces discretisation errors into physical quantities. 
These errors generally vanish with a positive power of 
the lattice spacing (a). 

Nonperturbative calculations in lattice QCD can be 
performed using Monte Carlo methods. Lattice arti- 
facts can be removed by reducing the lattice spacing 
used in numerical calculations. However, the com- 
putational cost increases as 1/a 7 (keeping the other 
parameters fixed). Alternatively, one can reduce dis- 
cretisation errors by adding higher-dimensional oper- 
ators to the action. This is called improvement. With 
improved actions the computational effort needed to 
perform reliable lattice QCD calculations can poten- 
tially be significantly reduced. This idea is behind 
much of the important progress made in lattice QCD 
calculations in recent years and has been an increasing 
part of research in lattice field theory. 

The main obstacle for obtaining quantitative results 
(at the few percent level) from numerical simulations 
of lattice QCD has always been the computational ef- 
fort associated with the proper inclusion of sea quark 
effects. Several years ago, substantial progress was 
made on this problem, in large part due to the de- 
velopment of an improved staggered fermion action 
@. For the first time, computationally efficient lat- 
tice simulations with realistic sea quark effects have 
become possible. 



1.2. Light Quark Methods 

The simplest lattice quark action replaces the co- 
variant derivative in the continuum action by a dis- 
crete difference operator. This so-called naive action 
suffers from the doubling problem. For every con- 
tinuum quark flavor, it has 15 additional unphysical 
flavors, called tastes. The staggered quark action com- 
bines four of these tastes into one Dirac field, by stag- 
gering the quark fields on a hypercube. This leaves 
four unphysical flavors (tastes). This action suffers 
from large 0(a 2 ) lattice spacing artifacts due to taste 
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changing interactions. The Asqtad action is an im- 
proved staggered action where all tree-level discretisa- 
tion errors are removed Q . Its leading lattice spacing 
errors are therefore of 0(a s a 2 ) and greatly reduced 
compared to the original staggered action. The Asq- 
tad action is the computationally most efficient light 
quark action available. However, in order to use it 
for the sea quarks in numerical simulations, the un- 
physical flavors must be removed. The sea quarks are 
present in the fermion determinant of the path inte- 
gral. To simulate two degenerate light (up and down) 
sea quarks, the MILC collaboration simply takes the 
square root of the light quark fermion determinant. 
For the strange sea quark, they take the fourth root 
of the determinant. This procedure is still somewhat 
controversial, but there is a growing body of evidence 
that its effects are controllable and disappear in the 
continuum limit 3]. The Asqtad action with the 
square root trick has been extensively tested in nu- 
merical simulations, most prominently in Ref. 0]. 

The HISQ (highly improved staggered quark) ac- 
tion is another version of an improved staggered action 
Q. Like the Asqtad action, it removes all tree-level 
0(a 2 ) errors. The 0(a s a 2 ) errors in the Asqtad action 
are rather large, due to taste changing interactions 
which appear at one-loop order. The HISQ action re- 
duces the 0(a s a 2 ) taste-changing effects by roughly a 
factor of three over the Asqtad action. The HISQ ac- 
tion has not yet been used to generate nf = 2+1 sea 
quark ensembles. Its computational cost is expected 
to be about a factor of two larger than the Asqtad 
action. 

Other light quark methods include the Wilson ac- 
tion and its improvements J6[ , Domain Wall Fermions 
and Overlap fermions [J with increasing compu- 
tational cost. The Wilson action solves the doubling 
problem by adding a dimension five operator which 
breaks chiral symmetry. Domain Wall Fermions solve 
the doubling problem by adding a fifth dimension, 
while keeping chiral symmetry almost exact. Overlap 
fermions have exact chiral symmetry, but a compli- 
cated operator structure. 

1 .3. Heavy Quark Methods 

On the lattice, heavy quarks with am,Q large, are 
best treated within an effective field theory framework 
(NRQCD or HQET). One can start with an effective 
field theory, and discretise it as in Ref. [ll|, for ex- 
ample. Alternatively, one can start with a relativistic 
lattice action and analyze its mass dependent discreti- 
sation errors using effective field theories. The charm 
quark it too light for a straightforward implementa- 
tion of the former approach, so we will focus on the 
latter. 

The Fermilab approach [12] starts with the im- 
proved relativistic Wilson action Q and the observa- 



tion that the Wilson action has the same heavy quark 
limit as QCD. With a simple prescription, the Wilson 
action can be used for heavy quarks without errors 
that grow with the heavy quark mass, (amg)". This 
approach can be used for both charm and and beauty 
quarks. With the improved Wilson action, the leading 
discretisation errors are 0(a s A/rriQ) and 0(A/itiq) 2 . 

The HISQ action is so highly improved that it can 
be used for charm quarks with an additional tuning 
of a parameter in the action, provided that the lat- 
tice spacing is small enough J5J. The leading mass 
dependent discretisation errors are formally of order 
0(a s (am c ) 2 ) and 0(am c ) 4 . 

1.4. Systematic Errors 

The most important sources of systematic error in 
lattice QCD calculations are sea quark effects; using 
unphysically large masses for the up and down quarks; 
discretisation effects; finite volume effects; and renor- 
malisation effects. 

In order to be phenomenologically relevant, a lattice 
QCD calculation must use gauge configurations that 
include the effects of three light sea quarks. Since 
the masses of the up and down quarks are generally 
taken to be degenerate, this is also referred to as the 
nf = 2+1 case. 

Until roughly five years ago, almost all lattice QCD 
calculations used ensembles generated either in the 
quenched approximation or with an incorrect number 
of sea quarks (generally, rif = 2) because of the com- 
putational cost associated with including sea quarks in 
the simulations. The quenched approximation omits 
sea quark effects entirely, at the cost of adding a sys- 
tematic error in the range of 10 — 30% for physical 
quantities involving stable hadrons 0]. This error 
must be determined on a case by case basis. Simu- 
lations with an incorrect number of sea quarks carry 
a similar systematic error, which is hard to estimate 
a priori. 

The computational cost increases with decreasing 
sea quark mass as m^ 2 ' 5 . All simulations to date use 
masses for the light sea quarks which are larger than 
the physical up and down quark masses. (Note, the 
strange quark mass is large enough to be simulated 
at its physical value.) We can use chiral perturba- 
tion theory (ChPT) to guide the extrapolations from 
the light sea quark masses used in the simulations to 
the physical masses. ChPT is an effective theory of 
QCD, which can be applied to (lattice QCD calcula- 
tions involving) pions and kaons. It can be combined 
with heavy quark effective theory and be applied to 
heavy-light systems, such as D and B mesons. Fur- 
thermore, it can be extended to include the leading 
light quark discretisation errors. Indeed, this has been 
done for the taste changing interactions of the Asqtad 
action and is called staggered ChPT (SChPT) |lo[. 
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The ChPT extrapolations are a significant but con- 
trollable source of systematic error. In order to keep 
this error at the few percent level or less, one needs to 
include simulations with a range of light sea masses, 
keeping mi < m s /2. The lattice QCD calculations 
described here include light sea quarks with masses in 
the range mi = 1/10 m s - 1/2 m s . Hence, the lightest 
light sea quark masses are only a factor of 2 — 3 away 
from their physical value. 

As described in the previous sections, the lattice 
actions give rise to discretisation errors. They can 
usually be estimated a priori using power counting 
arguments. However, even with improved actions, it 
is important to study and possibly remove these errors 
by repeating the calculation at several lattice spacings. 

1.5. Simulation Parameters 

The simulation parameters of the nj = 2 + 1 sea 
quark ensembles generated by the MILC collaboration 
using the Asqtad action (with the square root trick) 
are listed below and are shown in Figure [1] Each en- 
semble contains between 450-800 configurations. The 
ensembles contain one sea quark with a mass near the 
strange quark mass, m Sl and two degenerate light sea 
quarks with masses, mi. 



• a = 0.15 fm; 
0.6 m s . 



mi = 0.1m s , 0.2m s , 0.4m s 



• a = 0.12 fm; mi = 0.125 m s , 0.25 m s , 0.5 m s , 
0.75 m s . 

• a = 0.09 fm; m; = 0.1 m s , 0.2 m s , 0.4 m s . 




0.04 0.08 0.12 0.16 

lattice spacing, a (fm) 



Figure 1: Simulation parameters for the MILC ensembles 
with nj =2 + 1. showing mi/m s vs. lattice spacing a (red 
squares). The physical point is at mi/m s = 1/25 (pink 
burst). 



2. Semileptonic D Meson Form Factors 

The semileptonic decays D — > K{ii)lv are medi- 
ated by weak vector currents. The hadronic matrix 
elements for semileptonic decays are parameterized in 
terms of form factors. In our case there are two form 
factors, conventionally f+{q 2 ) and fo(q 2 ). The form 
factors are functions of the virtual W boson momen- 
tum transfer, q 2 , or, equivalently, the recoil momen- 
tum of the daughter meson. This introduces addi- 
tional lattice spacing errors: 

(AlV„|£>) tat = (K\V^\D) cont + 0(a PK ) n (2) 

Hence, discretisation errors are smallest, when pk is 
small and q 2 ~ <7m ax = {mr> — mx) 2 - 

The finite lattice volume provides an infrared cut- 
off, and therefore a minimum value for finite momen- 
tum, p m in = tt. Lattice three-momenta can be writ- 
ten in terms of p min as p = p m m(n X t % f n z ), where 
n x ,n y ,n z are integers. For example, for a = 0.1 fm, 
L = 20, p min = 620 MeV. 

To date, the only lattice results for semileptonic D 
meson form factors with rif = 2+1 are from the Fermi- 
lab Lattice and MILC collaborations [13j. They use 
the MILC a — 0.12 fm lattices with light sea quark 
masses in the range mi — l/8m s - 3/4m s , the Asqtad 
action for the light valence quarks and the Fermilab 
action for the charm quark. Staggered chiral pertur- 
bation theory is used to extrapolate to the physical 
light quark masses and to remove the leading discreti- 
sation errors due to taste violations. 

Figure [5] shows a comparison of the lattice QCD re- 
sult for the normalization f+(0) for D — > Klv with 
experimental determinations (where V cs is taken from 
other sources). The results are in very good agree- 
ment; however, the lattice QCD result has much larger 
errors than the experimental determinations. The 
comparison between lattice theory and experiment for 
fl(0) is similar Q. 

The shape of the form factor can also be determined 
in lattice QCD. However, in Ref. [l3j the form factors 
were calculated at a few values of recoil momentum. 
Then the BK parameterisation was used to deter- 
mine the q 2 dependence of the form factors. Since the 
errors increase with recoil momentum, the shape of 
the form factors is fixed mainly from the form factors 
near q^ ax and from using the BK parameterisation 
[l5j]. The lattice QCD result appeared before the new 
measurements by the FOCUS [l7[ and Belle [l~8] col- 
laborations were announced, so it is one of a very few 
lattice QCD predictions. Figure [3] [l5( | shows a com- 
parison of the lattice prediction for the q 2 dependence 
with experimental data from the Belle collaboration 
[l8j . The agreement is excellent. However, a quanti- 
tative comparison between the BK shape parameter 
determined from experiment and lattice theory is dif- 
ficult to interpret, as eloquently argued by Richard 
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Figure 4: The form factors for B — > 7vlv vs. q 2 from fit- 
ting to the z-expansion. The preliminary results are from 
Ref. fH. 




Figure 2: f+{0) in comparison from Ref. [l4| 




Hill 19]. A model independent parameterisation of 
the shape based on the z-expansion would avoid this 
difficulty 20]. The z-expansion is being used by the 
Fermilab Lattice and MILC (FNAL/MILC) collabo- 
ration to parameterize the q 2 dependence of the form 
factors for B — » irlv [21]. This works quite well, 
as shown in Figure [U Any new lattice analysis of 
semileptonic D meson decay form factors will (should) 
adopt the z-expansion to determine the shape. 

A number of additional improvements are possible 
in future calculations. Twisted boundary conditions 



can be used to adjust the lattice momenta to arbitrar- 
ily small values [22j, which would improve the shape 
determination. Double ratio methods similar to what 
has been developed for lattice studies of B — > Dlv | 23] 
can be adapted to semileptonic D meson decays [241 ] . 
This may lead to reduced statistical errors as well as 
improvements of some of the systematic errors. 



3. Leptonic Decay Constants f D and f Ds 

Charm leptonic decays provide another important 
test of lattice QCD. The lattice methods for calculat- 
ing decay constants in the charm and beauty meson 
systems are the same. Indeed, with the Fermilab ap- 
proach one uses the same heavy quark action in both 
systems and the heavy quark discretisation errors are 
expected to be larger for D mesons than for B mesons. 

There are now results from two groups 
(FNAL/MILC and HPQCD). Both use the MILC 
ensembles at a = 0.09 fm, 0.12 fm, 0.15 fm. 

The first FNAL/MILC results came out in 2005 
[25| , just days before CLEO-c announced its first pre- 
cise determination of fjj+ [261 ] ; the two results were in 
good agreement. 

The HPQCD collaboration announced their results 
for decay constants with much reduced errors this 
summer [27j and FNAL/MILC presented updated re- 
sults at the Lattice 2007 conference [2^], also with 
reduced errors. The new FNAL/MILC analysis was 
done "blind" , where an overall unknown offset was 
added to the lattice data. The final results were un- 
blinded shortly before they were presented at the Lat- 
tice 2007 conference, making this the first (intention- 
ally) blind lattice analysis. Table [I] compares the main 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



5 



features of the two calculations. More details about 
the HPQCD and FNAL/MILC calculations, includ- 
ing discussion of the error analysis and plots of chi- 
ral and lattice spacing extrapolations can be found 
in Refs. [H and (H respectively. The FNAL/MILC 
analysis includes more lattice ensembles, more valence 
quark masses per ensemble, and uses staggered chiral 
perturbation theory (Staggered ChPT) to remove the 
leading light quark discretisation errors. The HPQCD 
collaboration considers only the case m q = mi , where 
m q denotes the light valence quark mass and m; de- 
notes the light sea quark mass. They use continuum 
ChPT with generic 0(a 2 ) terms added in the chiral 
fits. 

The main difference between the two calculations 
is the valence quark actions. The HPQCD collabo- 
ration uses the HISQ action for all (charm, strange 
and light) valence quarks, whereas the FNAL/MILC 
collaboration uses the Fermilab action for the charm 
quarks and the Asqtad action for the strange and light 
valence quarks. Since the HISQ action is more im- 
proved than the Fermilab action, the HPQCD result 
has much smaller heavy quark discretisation errors. 
This is the main reason for the difference in the total 
errors between the two results. 

Table |H] compares the error budgets for the 2005 
FNAL/MILC calculation with the FNAL/MILC Lat- 
tice 2007 one. The error reduction is mainly due to 
including three MILC ensembles at a = 0.09 fm (and 
8-12 different valence masses). This reduces the heavy 
quark and light quark discretisation errors, and better 
constrains the staggered ChPT. 

Figures [6j and [7] compare the lattice results for 
Jd+i fi) B i an d /d s //d+, respectively, to the corre- 
sponding experimental averages. The experimental 
averages are from Ref. [14]. The new CLEO-c result 
fo s = 275 ± 10 ± 5 presented at this conference by 
Steven Blusk [3(|, is very similar to Ref. [Til ]. 

The FNAL/MILC results agree with the experimen- 
tal averages at the one sigma level. The HPQCD re- 
sults agree very well with the FNAL/MILC results. 
There is a hint of disagreement between the HPQCD 
result for fjj 3 and the experimental average at the 
two sigma level. However, the experimental determi- 
nations of the decay constants must assume a value 
for the CKM angle V cs from other sources. We are 
approaching a level of precision, where tests of lat- 
tice QCD should be performed on CKM free quanti- 
ties such as the ratio of semileptonic to leptonic decay 
rates suggested in Ref. [31| . 



4. Conclusions and Outlook 

With the generation of the MILC ensembles, the 
stakes for lattice QCD calculations have risen. We are 
now able to calculate the simplest quantities to a few 
percent accuracy. As always, repetition is desirable 
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Figure 5: Comparison of lattice QCD results for f D + with 
experiment. 
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Figure 6: Comparison of lattice QCD results for fr>, with 
experiment. 



to test different lattice methods against each other. 
To date, all lattice calculations that include realistic 
sea quark effects use the MILC ensembles with rooted 
Asqtad sea quarks. As mentioned in section [l~2l the 
Asqtad action carries the smallest computational cost 
of any light quark action. Nevertheless, recently other 
collaborations have started to generate ensembles with 
different sea quark actions. An overview is given in 
Figure [H It shows that the other ensembles are be- 
ing generated at similar values of lattice spacing and 
light quark masses as the MILC ensembles. The MILC 
collaboration continues to generate new ensembles at 
even smaller lattice spacings. They are also generating 
additional configurations for the existing ensembles to 
further reduce statistical errors. As in experiment, in 
lattice QCD smaller statistical errors give better con- 
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Table I Comparison of the main features of the HPQCD and Fermilab Lattice/MILC calculations. 
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HISQ action for charm and light valence 


Asqtad action for strange and light valence 
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a (fm) 


mi/m 3 sea quark 


a (fm) 


mi/m s sea quark 


0.09 
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1 valence quark mass/ensemble, m va i e nce = m sea . 


Partial nonperturbative renormalisation 


Nonperturbative renormalisation from PCAC 


Staggered ChPT fits to all valence 


Continuum ChPT +0{a 2 ) terms fit to all 


and sea quark ensembles together 


ensembles together 


Blind analysis for Lattice 2007 







Table II Comparison of the error budget of the 2005 
FNAL/MILC results with the Lattice 2007 results. Num- 
bers are given in percent. 
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0.5 


2.7 
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2.8 
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3.4 
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higher order PT 
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0.3 




+ other small sources (finite volume, . . 


) 
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trol over systematic errors. Hence, lattice calculations 
based on the MILC ensembles will continue to become 
more accurate. 

Any modern lattice QCD calculation that claims 
phenomenological relevance must include a serious 
systematic error analysis. To be relevant, it must in- 
clude the correct number of sea quarks (which all the 
ensembles shown in Figure [5] do) . While the masses 
of the light sea quarks are still unphysically large, 
it must also include a study of the light quark mass 
dependence with sufficiently small sea quark masses. 
Among other sources of error, discretisation effects 
must be estimated and tested by repeating the cal- 
culation at more than one lattice spacing. 

In summary, we should expect to see lattice results 
from these new ensembles in the near future. They 
will provide important consistency tests of the lat- 
tice methods, and in particular of the square root 
trick used by the MILC collaboration to generate their 
rtf = 2+1 ensembles. 



Figure 7: Comparison of lattice QCD results for fo B / fo+ 
with experiment. 
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